Charge excitations in LiV205 and NaV205: similarities and differences 
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We calculate the optical conductivity of LiV2 05 and NaV2 05 using exact numerical diagonal- 
ization of a quarter-filled extended Hubbard model on a system of coupled ladders. In particular, 
electronic correlations are treated exactly, and a quantitative agreement between calculated and ex- 
perimentally observed optical conductivity of these two vanadium oxides is presented. Furthermore, 
it is found that LiV205 differs from NaV205 not only in the charge ordering pattern but also in 
the nature of the inter-ladder coupling: In contrast to LiV205, in NaV205 neighboring ladders are 
coupled by a strong Coulomb repulsion, and not by inter-ladder hopping. 
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In the past years, low-dimensional transition metal 
compounds have been intensively investigated because 
of their unconventional spin and charge excitation spec- 
tra In this respect, the vanadate Q;'-NaV205 has 
attracted particular interest as a ladder system at quar- 
ter filling, containing only one equivalent V site with a 
formal valence -1-4.5 |^,^. The magnetic susceptibility 
of NaV205 can be well described by a S" = 1/2 antifer- 
romagnetic Heisenberg chain with exchange interactions 
of J = 440 and 560 K for temperatures below Q and 
above |^ the transition temperature Tc w 34 K. The 
low-temperature phase is found to be charge ordered , 
but the nature of this transition is still under discus- 
sion 0-^. The much less studied 7-LiV205 belongs to 
the same family of vanadium oxides and exhibits a one- 
dimensional S = 1/2 Heisenberg like behavior with an 
exchange interaction of J = 308 K |p^ , |rT| . In contrast 
to NaV205, there is no indication of a phase transition 
at lower temperature. Since both compounds are struc- 
turally related it is interesting to clarify the microscopic 
origin of the different physical properties of both vana- 
dates. This will be discussed in this paper on the basis 
of the optical conductivity. 

Recently, the optical conductivity of LiV2 05 and 
NaV205 has been measured ||l2|. In the energy range 
from to 3 eV similar peaks in the E || a spectra of both 
materials were found. On the other hand, a complete 
suppression of the peaks in the E || 6 spectrum of LiV2 05 
was observed. In this paper, we show that this suppres- 
sion results not only from the double chain charge order- 
ing pattern but also from the strong inter-ladder hopping 
in LiV205. Consequently, it is found that both vanadates 
differ not only in the charge ordering pattern but also 
in the nature of the inter-ladder coupling. Recently, we 
studied the optical conductivity of NaV205 so that 
in the present paper we concentrate on LiV205 and on 
the comparison of both vanadium oxides. 

At room temperature LiV205 and NaV205 have or- 
thorhombic crystal structures that are described by space 
groups Pnma and Pmmn, respectively |14,15[. Both 



compounds consist of layers of VO5 square pyramids (see 
Fig. 0). In contrast to NaV205, a structural analysis 
p4| for LiV2 05 shows two inequivalent vanadium sites 



which were also found by NMR experiments |16|. These 
two vanadium sites were assigned a valence of V^"'' and 
V^+, respectively, and they form two different zig-zag 
chains along the b direction jl^. Therefore, magnetic 
V''^"'" double chains are separated by nonmagnetic 
double chains. Note that this charge ordering pattern 
is not comparable with the in-line order discussed for 
NaV205 [0. As can be seen from Fig. |l|, the VO layers 
are more corrugated in LiV205 than in NaV205 because 
the Li atoms are smaller than the Na atoms. 
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FIG. 1. Representation of the (a) NaV205 and (b) LiV205 
room temperature crystal structure, (c) shows a schematic 
projection of the LiV205 structure onto a {a,b) plane where 
only the vanadium ions are drawn. 

In Ref. [|l3| we studied the optical conductivity of 
NaV205 using a quarter-filled extended Hubbard model 
of the V 3(i electrons on the vanadium lattice which forms 
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a coupled ladder system. In the present study we modify 
this model for LiV205 by an additional on-site energy Si 

H = - ^ tij (4^Cj,^ + H.c.^ + "iT'^a 

+ ^ + (1) 

to cause the double chain charge ordering observed in 
LiV205 llj]. In Eq. denotes summation over 

all pairs of nearest neighbors, and spin a =T,i- cJcr 
are electron creation operators at vanadium sites, rii ~ 



is the occupation-number operator, and U 



denotes the Coulomb repulsion between electrons on the 
same site. The hopping parameters tij of Eq. ^ are de- 
fined in panel (c) of Fig. |l|. The inter-site Coulomb in- 
teractions Vij between nearest vanadium neighbors on a 
rung, on a leg or on different ladders are Va, Vb and Vxy 
The on-site energy is -A/2 (A/2) for the V'^+ {Y'^+ 
sites [see panel (c) in Fig. In contrast to NaV205, in 
LiV205 the ladders are less strictly situated in an (a, b) 
plane (see Fig. |l]). Thus, the system shown in panel (c) 
of Fig. I has more 3D character than NaV2 05. 

For an insulating system, the optical conductivity 
CTq,([j) is given by the Kubo formula for finite frequency 
response 



a„(w) = -Re / dte*"*(0|j„(t)j„|0). 
Jo 



(2) 



Here ja with a = a,b are the components of the cur- 
rent operator parallel to a or & direction. Equation (|^) 
is valid for zero temperature, whereas the experiments 
have been carried out at finite temperatures |12[| . How- 
ever, it is known from both experiment and theory that 
at least the spectra for NaV205 depend only weakly on 
temperature [18|,|9|. Therefore, we may restrict ourselves 
to calculations for zero temperature. 

In the following we evaluate Eq. (||) using the standard 
Lanczos algorithm ||2^ which is limited to small clusters. 
Therefore, first we have to check if our results are suf- 
ficiently converged with respect to system size. Since 
the hopping between V"'''" ions in LiV2 05 is very small 
this material can be described by the double lad- 
der system shown in panel (c) of Fig. 0. Consequently, 
the system is infinite only in 6-direction, and we shall 
always use open boundary conditions for this direction 
(to enlarge the effective cluster size). However, one has 
to make sure that electrons on the edges of the cluster 
are still embedded in the local Coulomb potential that 
results from the double chain charge ordering. For this 
purpose, sites on the edges of the clusters that are oc- 
cupied in the perfectly ordered state need an additional 
on-site energy Vb to simulate the influence of occupied 
sites outside of the cluster. To prove that the results arc 
sufficiently converged with respect to system size we have 



evaluated the double ladder system shown in panel (c) of 
Fig. 1^ both with open and periodic boundary conditions 
and found only small differences (not shown). 

The distances within the vanadium ladder structures 
in LiV205 and NaV205 are similar. Therefore, we start 
from the parameters of the ladder used in Ref. jl^] 
for NaVaOs (ta - 0.38 eV, h = 0.17 eV, tdiag = 0, 
Va = 0.8 eV, Vb = 0.6 eV, U = 2.8 eV) to describe the 
optical conductivity of LiV205. Consequently, we have to 
determine only three parameter values. A, Vxy, and t^y 
For LiV2 05 the choice Vxy = is possible because the 
Coulomb repulsion between nearest neighbors of neigh- 
boring ladders counteracts the double chain charge or- 
dering caused by A. This means that the parameters A 
and Vxy are not independent from each other. 
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FIG. 2. Optical conductivity [panel (a): (Ta(cj), panel (b): 
at{LL>)] of a single ladder {t^y ~ 0) where the hopping param- 
eters from Ref. ^ (dotted lines, ta = 0.38 eV, U = 0.17 eV, 
idiag ~ 0) and from Ref. (dashed lines, ta ~ 0.35 eV, 

tt = 0.02 eV, tdiag ~ 0.1 eV) were used. In addition, the re- 
sults for the double ladder system shown in panel (c) of Fig. ^ 



{ta = 0.35 eV, tb = 0.02 eV, tdia 



0.1 eV, 



-0.18 eV) 



is plotted with full lines. The parameters of the Coulomb in- 
teraction are U = 2.8 eV, 14 = 0.8 eV, Vb = 0.6 eV, 14;/ = 0, 
and A — 1.7 eV. The theoretical line spectra are broadened 
with Gaussian function of width 0.1 eV. 



First we discuss the case of a single ladder {txy = 0) for 
which the optical conductivity is shown as dotted lines 
in Fig. ||. As can be seen from panel (a) of Fig. |^, (Tq is 
dominated by an excitation at 1.0 cV which corresponds 
to a transition from a bonding to an anti-bonding state 
of a singly occupied rung that has been discussed be- 
fore for N aV205 p^ . If one considers a single rung one 
obtains y^S'^ + for the excitation energy of this tran- 
sition where an effective on-site energy (5 = A — 2Vf, has 
to be used to take the chain charge ordering into ac- 
count. Therefore, the on-site energy A determines the 
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excitation energy (of the transition from a bonding to 
an anti-bonding state). A value of A = 1.7 eV can be 
obtained directly from the experimentally observed peak 
position [T^ . In contrast to a a, for we find two ba- 
sic mechanisms for excitations. The structures at 1.2 eV 
and 1.8 eV [see dotted line in panel (b) of Fig. ^ can 
be interpreted as a transition to states with one unoccu- 
pied and one doubly occupied rung [excitation energies 

(A + Vq - SVb) and (A + K - Vb)]. The excitations 
around 2.8 eV result from the creation of one doubly oc- 
cupied site (excitation energy ~ U). 

Next we discuss the influence of the inter-laddcr hop- 
ping txy on the optical conductivity. A moderate value of 
txy affects a a only weakly and leads, on the other hand, 
to a suppression of low-energy peaks in (Th (not shown) . 
This suppression becomes stronger if A increases. (Be- 
cause of different charge ordering pattern this influence 
of A is not identical with the findings in Ref. |2^.) One 
would obtain good agreement between the calculated and 
the experimentally observed p2| optical conductivity for 
txy = 0.3 . . . 0.4 cV. However, according to Ref. the 
magnetic properties of the model along 6-direction can 
be described by a single Hcisenbcrg chain along a zig-zag 
line. The choise of txy ~ 0.3 . . . 0.4 eV for the inter- 
ladder hopping supposed above would lead to a strong 
frustration |24| 



A 

Jxy 
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which is not in accordance with NMR [|Tl| and suscep- 
tibility measurements [0. However, this discrepancy is 
caused by neglecting the diagonal hopping tdiag^ Tight- 
binding fits to the band structure found in LDA cal- 
culations with nonzero tdiag lead to a small value of th 
both for NaV2 05 and for LiV2 05 [||,|l|. Therefore, we 
use in the following the set of tight-binding parameters 
from Ref. (2l) for LiVaOg {ta = 0.35 eV, 4 = 0.02 eV, 
^diag = 0.1 eV, txy = —0.18 eV). Since the diagonal hop- 
ping idiag is included, the frustration parameter (^) is 
now negligibly small [^. The Coulomb interactions U, 
Va, and Vb of the ladder arc not changed. In contrast 
to LiV205, for NaV205 we use a parameter set with 
^diag = as in Ref. |l3j to avoid strong finite-size effects 

From a comparison of the optical conductivity of a 
single ladder for both parameter sets (dotted and dashed 
lines in Fig. ^ one can conlude that the results are not 
fundamentally changed by the tight-bindung parameters 
from Ref. [|l|- But these new parameters lead to a signif- 
icantly decreased spectral weight of the low-energy struc- 
ture in (Tb because the mobility of a created doubly occu- 
pied rung is now more restricted due to the small value 
of tb. Nevertheless, the above analysis of the basic mech- 
anisms for excitations on a single ladder remains valid 
even for tdiag 7^ 0. 



Now we return to the influence of the inter-ladder hop- 
ping on the optical conductivity. For this purpose, we 
compare in Fig. |^ the results of a single ladder (dashed 
lines) with the optical conductivity for the complete dou- 
ble ladder system shown in panel (c) of Fig. |l] (full lines) 
where we use the hopping parameters from Ref. [ pl[ . As 
can be seen from panel (a) the excitation resulting from 
the transition between a bonding and an anti-bonding 
state of a rung discussed above depends only weakly on 
txy In contrast to (Ta, ab is changed strongly by a finite 
txy In particular, some new structures arise, and the 
spectral weight is divided among different excitations. 
Note that the influence of txy does not depend on its 
sign. 
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FIG. 3. Comparison of experimental data for LiV205 
[panel (b)] and for NaV205 [panel (d)], taken from Ref. |l^, 
and the results of the exact diagonalization [panels (a) and 
(c)] whereby aa (crt) is plotted with full (dashed) lines. The 
hopping parameters are taken from Ref. |l| (LiV2 05) and 
Ref. 1^, respectively. The Coulomb interactions on the lad- 
ders {U = 2.8 eV, Va = 0.8 eV, Vt = 0.6 eV) are the same 
for both compounds. Additional parameters are A = 1.7 eV 



for LiV2 05 and 14, 



0.9 eV for NaV205. The theoretical 



line spectra have been convoluted with a Gaussian function 
of width 0.1 eV. 

In Fig. ^ the calculated optical conductivities for 
LiV205 and NaV205 [panels (a) and (c)] are compared 
to the experimental spectra [panels (b) and (d)] from 
Ref. [|l^. We restrict the following discussion of the op- 
tical conductivity to the energy range up to 2 cV since 
transitions with a higher excitation energy involve oxygen 
orbitals [|l^. Good agreement for the optical conductiv- 
ity of NaV205 is found [compare panels (c) and (d) of 
Fig. ^ where the hopping parameters {ta ~ 0.38 eV, 
tb = 0.17 eV, txy = 0.012eV) and the on-site Hubbard 
interaction {U ~ 2.8 eV) are taken from Ref. |]. The 
values of the intersite Coulomb interactions Va ~ 0.8 cV, 
Vb — 0.6 eV and Vxy = 0.9 eV have been adjusted to 
obtain correct peak positions (for details see Ref. [Of). 
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Note that the loss function in EELS experiments for 
NaV205 can also be described with the same Hamilto- 
nian, and the same set of model parameters For 
LiV205 we obtain good agreement of the calculated op- 
tical conductivity [panel (a) of Fig. ^ and the exper- 
imental data of Ref. [panel (b) of Fig. ^. Again 
we use the hopping parameters from LDA calculations 
|l| {ta = 0.35 eV, h = 0.02 eV, tdiag = 0.1 eV, 
txy = —0.18 eV) and the Coulomb interactions of the 
ladder [U = 2.8 eV, K = 0.8 eV, Vb = 0.6 eV) found 
for NaV205. The on-site energy A = 1.7 eV is adjusted 
to obtain the correct peak position of a a- Thus wc use 
only one free parameter for LiV205. The effective on- 
site energy 5 = A — 2Vb = 0.5 eV that determines the 
charge ordering along the chains agrees well with the 
value 5 = 0.3 eV found in Ref. [|21||. To summarize 
the comparison of the calculated and the experimentally 
observed optical conductivity for LiV2 05 and NaV2 05 
(Fig. H) one can state that a good agreement is found for 
both vanadates. In particular, the complete suppression 
of the low energy structure in ai, of LiV205, and the peak 
shapes arc reproduced. 

Finally, wc want to compare the results for LiV2 05 
and NaV2 05 in some detail. As shown above, the lad- 
der parameters are similar for both compounds. This 
means that the electronic properties of the single lad- 
ders are similar in LiV205 and NaV205 even though one 
finds a double chain charge ordering in LiV205 caused 
by an on-site energy A. Apart from that, the main dif- 
ference between LiV205 and NaV2 05 is the character of 
the inter-ladder coupling. In LiV205 one finds a strong 
intcr-laddcr hopping t^y [keep in mind that this coupling 
docs not go beyond the double ladder systems shown in 
panel (c) of Fig. |^- which is very small in NaV2 05. On 
the other hand, in NaV205 the ladders are coupled due 
to a strong inter-ladder Coulomb interaction Vxy Note 
again that the choise Vxy = is reasonable for LiV2 05 
since A and Vxy are not independent from each other. 
The character of the inter-ladder coupling determines the 
low energy features of Ub'. the well defined structure in 
NaV205 is suppressed in LiV205 (see Fig. ||) where a 
finite txy leads to new excitations among which the spec- 
tral weight is divided. 

In conclusion, we have obtained a quantitative de- 
scription of the highly anisotropic optical conductivity of 
LiV205 and NaV205 by exact diagonalization of small 
clusters. In particular, the complete suppression of the 
low-energy peaks in the spectrum of LiV2 05 with elec- 
trical field parallel to the leg direction is reproduced. 
The good agreement of calculated and experimental op- 
tical conductivity for both vanadates implies that the 
electronic properties of the single ladders arc similar in 
LiV205 and NaV205 (even though one finds a double 
chain charge ordering in LiV2 05). The main differences 
between the optical properties of LiV2 05 and NaV2 05 
are due to the completely different nature of the inter- 



ladder coupling. In contrast to LiV2 05, in NaV2 05 dif- 
ferent ladders arc coupled by a strong Coulomb repulsion 
Vxy and not by inter-ladder hopping txy Because of the 
on-site energy A LiV205 can be viewed as an asymmetric 
quarter-filled ladder compound. Therefore, it would be 
interesting to compare with other experiments probing 
charge excitations like EELS or XPS. 
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